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Abstract — The objective of this paper is to design optimal sig- 
nature matrices for binary inputs. For the determination of such 
optimal codes, we need certain measures as objective functions. 
The sum-channel capacity and Bit Error Rate (BER) measures 
are typical methods for the evaluation of signature matrices. 
In this paper, in addition to these measures, we use distance 
criteria to evaluate the optimality of signature matrices. The 
Genetic Algorithm (GA) and Particle Swarm Optimization (PSO) 
are used to search the optimum signature matrices based on 
these three measures (Sum channel capacity, BER and Distance). 
Since the GA and PSO algorithms become computationally 
expensive for large signature matrices, we propose suboptimal 
large signature matrices that can be derived from small sub- 
optimal matrices. 



I. Introduction 

CODE Division Multiple Access (CDMA) is a method for 
reliable data communication among multiple users and 
is the standard 3G wireless systems. The general model of a 
CDMA system is defined as: 



Y ^ AX + N 



(1) 



where A is the m x n signature matrix, m is the length of 
signatures and n is the number of users. X is an n x 1 user 
column vector and N is the Additive White Gaussian Noise 
(AWGN) vector N =[Ni,..., iV„j]^, such that A^/s are i.i.d 
random variables. For binary input CDMA system, entries of 
X are binary {i.e., {±1}), with uniform distribution. 

Due to the bandwidth constraint of mobile communication 
systems it is desirable to use more user than possible. Thus, 
an overloaded CDMA comes to life when the number of users 
exceeds the length of the signatures(r7i > in this situation 
orthogonal signature matrices, such as Hadamard codes, can 
no longer be used. 

An loading factor for such systems is defined as follows: 



m 



(2) 



Most of the work in the evaluation of sum capacity has been 
done for large scale CDMA systems (asymptotic results) ||T|- 
|j3j|. For finite scale systems the actual capacity is not known, 
however there are lower and upper bounds for the sum capacity 
Q- ||6). A review of these papers is done by Hosseini in ||7). 

Pad et al. |4| presented optimum signature sets for binary 
input and CDMA matrices such as Codes for Overloaded 
Wireless (COW) matrices. Furthermore, a new ML decoder is 
introduced for large scale COW matrices. Alishahi et al. |5| 



evaluated the upper and lower bounds of sum channel capacity 
for binary CDMA systems with and without noise. Alishahi et 
al. |j6| generalized these previous works Q and p) to finite 
none binary descrete input and matrix entries. InPSl we got 
some partial results for optimal real matrices base on GA. 

In the present paper we aim to derive optimum real or binary 
signature matrices with binary input. For the optimality we 
use different measures such as sum channel capacity, bit error 
rate BER and in addition to them we use a distance criteria to 
reduce the computational complexity. It is noted that we can 
use these criteria for non binary input CDMA systems. 

Our paper includes the following main contributions: Firstly, 
we propose three main criteria for optimizing signature ma- 
trices; the channel capacity maximization, BER minimization 
and distance criteria. Secondly, Since it is difficult to evaluate 
these three criteria for large scale systems, we present a 
method to derive sub-optimum signature matrices by enlarging 
low dimensional ones. 

The rest of this paper is organized as follows: Optimiza- 
tion measures and methods for optimum signature matrices 
derivation are discussed in Section |ll] In Section |lll] numerical 
and simulation results based on GA and PSO are presented 
and compared Section |V] proposes a method to derive sub- 
optimum high dimensional signature matrices. Finally, Section 
VI concludes the paper and highlights the future works. 



II. Signature Matrix Optimization Measures 

In this section we will discuss three measures for designing 
sub optimal signature matrices. We will show that there is a 
trade of between accuracy and computational complexity of 
these criteria. Each criterion has certain properties which is 
useful in specific conditions such as low SNR values and high 
overloading factors. 

A. Channel Capacity Criterion 

The most precise method to derive optimized signature 
matrices is by using the sum capacity criterion. The sum 
capacity for a given dimension nxm and noise with variance 
(Ttv can be defined as: 

C{n,m,(Tis[) = max C(n, m, crjvl A) (3) 

where C{n, m, (Jn\A) is the sum capacity for a specific matrix 
A. Considering that A is deterministic in ^ and (j3]l, we have: 

C(n,m,aN\A) = max I(X:Y) = maxh(Y) - h(N) (4) 

P{X) P{X) 



2 



where h{y) is the differential entropy. According to the con- 
jecture mentioned in f9l and [5|, h{y) is maximized when X 
is uniform. For input uniform distribution since iV^'s are i.i.d. 
and thus f^iN) = YiiLi fni{ni), the probabihty distributed 
function (pdf) of Y is as bellow: 



E 

xe{±i}"> 




(5) 



where i/i and are the i entry of Y and the i row of 
matrix A, respectively. As a result of this, h{Y) is computable 
from (|5]l as: 



hriY) 



/y(y)log2(/r(y))dyi...d2/„ 



Consequently C(n, m, cttvI A) is derived from the mutual 
entropy in Q. The per-user sum capacity parameter is used 
instead of channel capacity, which can be defined as: 



(6) 



^Channel 

C = 

n 

where c is the normalized of C. It must be mentioned that the 
computational complexity of capacity measure is of 0(to") 
which results in an NP-hard algorithm. 

B. BER Criteria 

To compute the BER measure, a large array of bits (10^ 
bits) are produced, encoded, transmitted through the simulated 
channel and decoded at the receiver by an ML decoder. After 
the decoder, the probability of error statistically measured. 
We conjecture that the lower the BER of a signature matrix 
is, the higher its channel capacity will be; this conjecture is 
verified with simulation results. Hence, optimum matrices can 
be derived by minimizing the BER measure. 

The BER requires much less computation than the capacity 
evaluation. However both criteria are not practical due to ex- 
tensive computation. Therefore, we propose distance measures 
as discussed bellow. 

C. Distance criterion 

In this subsection, we propose three different methods based 
on the output constellation points (Zi), i.e.. The output points 
which are defined in the absence of noise: 



(7) 



Fig. [T] shows an example for constellation of a 2 x 3 real 
signature matrix (which is optimized using channel capac- 
ity method for E^/Nq = 5dB). This constellation is the 
projection of ±1 three dimensional cube on a plane. Since 
X,, -X, G {±l}"^i for every i = 1, 2, • • • 2" we have a sim- 
ilar symmetry property for Z/s. Therefore, the computational 
complexity is reduced by considering only half of Z^'s in the 
distance criterion. One way for optimization is maximizing 



Fig. 1 . Constellation of a 2 X 3 signature matrix optimized for Ei, / Nq = 5dB 
using the channel capacity criteria 



the minimum distance of the output constellation points. This 
method guarantees an upper bound for the probability of error 
for high values of SNR. The following equation explicitly 
defines Minimum Distance (MD) criterion: 



MD 



mm 



\Z, 



where \\U\\ represents the Euclidean norm of the vector 
U. Note that MD must be maximized in order to find a 
sub-optimum matrix. By taking a more analytical approach, 
suppose that xi's (the i*'^ element of the input vector X E 
{±1}"^^) are independent, hence the error probability of a 
block (with size of m) is derived as follows: 



2" 

p, = -yp 

2" ^ 

1=1 



u 

fc=l 



\Y-Z,\\' > WY-Zk 



(8) 



where Yi — Zi+N and Zi is the output vectors in the noiseless 
channel. Since the noise vector elements are i.i.d. Gaussian 
random variables with variance ctat, the upper bound for error 
probability can be calculated as: 



1 



2" 2" 



\Zi — Zj 
2a N 



(9) 



where Q{x) refers to Cumulative Distribution Function (CDF) 
of Gaussian distribution. Instead of maximizing MD, we 
minimize the upper bound given in (|9| which is equivalent 
to minimizing the following function: 



QD 



2" 2" 

EE 

1=1 j=i 



Q 



\Zi — Zj 
2ct 



(10) 



The minimization of the upper bound can give the better 
results than MD; however it is more computationally intensive. 
To reduce the computational complexity of Q{x), the follow- 
ing approximation can be used the following approximation: 



Q{x) 



0.7 exp - 




(11) 
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Thus we can use the following distance denoted by ED: 



2" 2" 



ED= l^l^ exp 



i=i j=i 



1.6 



(12) 



For high SNR values, ( 12 1 is simplified to a single exponential 



element which presents a similar behaviour for MD. The next 
section is a brief description for GA and PSO. 



III. Optimization Techniques 

In this section, we will describe the optimization techniques 
to find good signature matrices based on the criteria discussed 
on Section [n] that are applied to optimize the criteria. We use 
two known optimization techniques, namely GA and PSO. We 
show that since GA and PSO are methods for minimizing an 
arbitrary cost function, treating the criteria as cost functions, 
will lead to derive optimum signature matrices. 



A. Genetic Algorithm 



The Genetic Algorithm (GA) 1 12| employs the principal of 
survival of the fittest in its search process to select and generate 
individuals (design solutions) that are adapted to their environ- 
ment (design objectives/constraints). Therefore, over a number 
of generations (iterations), desirable traits (design characteris- 
tics) will evolve and remain in the genome composition of the 
population over traits with weaker undesirable characteristics. 
The GA is well suited and has been extensively applied to 
solve complex design optimization problems because it can 
handle both discrete and continuous variables with nonlinear 
objective and constraint functions. In this work we apply the 
GA to find sub optimum signature matrices based on criteria 
of sum capacity, BER, and the distance. The next algorithm 
is an alternative algorithm that can converge to some optimal 
matrices with the faster rate of convergence. The parameters 
and options of the GA are presented in Table U 



TABLE I 

Simulation parameters set for GA 



Population 



Size 
Type 

Creation function 

Lower bound 
Upper bound 

New Generation 



20 individual matrices 
Double vector 

Unifomi in first run Best results of 

previous runs afterwards 

-I 

+1 



Elite count 2 

Crossover fraction 0.8 

Migration direction Forward 

Mitigation factor 0.2 

Stopping criteria 



Number of iterations 
Function tolerance 



100 

10-6 



B. Particle Swarm Optimization 

PSO |[13|, |[I4) similar to GA is a computational method that 
optimizes a problem by iteratively trying to improve a candi- 
date solution, which results in an objective function. In PSO, a 
set of randomly generated solutions (initial swarm) propagates 
in the design space towards the optimal solution over a number 
of iterations (moves) based on large amount of information 
about the design space that is assimilated and shared by all 
members of the swarm. The PSO algorithm considers some 
candidate solutions in the search domain. During each itera- 
tion, the cost function of each candidate solution is calculated. 
Each candidate solution can be considered as a particle moving 
toward the minimum value of the cost function. As the first 
step, PSO chooses the candidate solutions randomly inside the 
search space. It should be mentioned that the PSO does not 
have any prior information about the cost function; it does not 
know which particles are near or far from the global minimum 
of the cost function. What PSO does, is to evaluate the cost 
value of each particle and just work with the corresponding 
cost values. The position of a particle is composed of its 
candidate solution, cost and velocity. Moreover, it remembers 
the least cost (the best fitness) that it has had thus far during the 
operation of the algorithm, called the individual best fitness. 
The candidate solution corresponding to this fitness is referred 
to as the individual best candidate solution or the individual 
best position. At last, the PSO seeks for and finds the least cost 
among all the particles in the swarm, named the global best 
fitness. Bellow we will compare the GA and PSO in terms of 
convergence rate. The simulation parameters and options set 
for the PSO algorithm are listed in Table |ll] 

TABLE II 

Simulation parameter set for PSO 



Particles position 

Number of particles 
Initial position 
Lower and upper bounds 
Particle Velocity 

Stopping criteria 



20 individual vectors 
Random, uniform distribution 
[1, +1] 

initiated randomly with uniform 
distribution 



Number of iterations 



100 



C. Convergence evaluation of GA and PSO 

We check the convergence behavior of the discussed opti- 
mization algorithms before we apply them for our analysis. 
It is important for GA to converge to the minimum point for 
each method. We can show the convergence rate by comparing, 
the best and the mean fitness values for every iteration that 
are the minimum and the mean of objective function in each 
iteration, respectively. For instance, considering MD as an 
objective function for n = 5 and rn = 4, there is no difference 
between the best and the mean values after 50 iterations (as 
shown in Fig. [2]), therefore, one concludes convergence of 
the GA . Furthermore, the variation of the mean fitnesses in 
various iterations shows that the GA explores almost the whole 
of the feasible populations, which implies global minimum 
as oppose local minimum. The best and the mean fitness 
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Fig. 2. The best value and the mean value of populations versus GA iterations 
for the MD method (/3 = |). 



values for each iterations for the PSO algorithm is depicted in 
Fig. [3] using MD as the objective function when n = 5 and 
TO ~ 4. In this figure unlike Fig. |2] the beast and the mean 
do not cross each other, however the best value is unchanged 
after a certain iteration. This figure shows that the algorithm 
converges after 28 iterations. A comparison between Fig. |2] 
and Fig. [3] demonstrates that the PSO converges faster than 
the GA. 




Fig. 3. The best value and the mean value of populations versus PSO 
iterations for the MD method (fi = |). 



IV. Numerical and Simulation Results 

This section presents the numerical and simulation results of 
the GA and PSO for proposed criteria, and compare GA and 
PSO in term of optimality. Sub-optimum signature matrices 
based on various criteria using GA and PSO are presented in 
the Appendix. 



A. Results of GA for Real Valued Signature Matrices 

In this subsection, we apply the GA to verify our results 
discussed in Section |ll] namely BER, channel capacity and 
distance methods. Firstly, we will compare distance methods 
with each other and then, we will show the simulation and 
numerical results of various methods. 



1 ) Comparison of Distance Criteria: In order to compare 
distance criteria in term of capacity. We obtain sub-optimal 
matrices for an arbitrary Ef, /N^ value and compare the capac- 
ity of these matrices for that E^/Nq value. Fig. |4] illustrates 
the per-user capacity of optimized matrices using distance 
criteria for different E^/No values. The matrices exhibit a 




Fig. 4. Per-user capacity of sub-optimum matrices based on distance methods 
versus different Ef,/No values (/3 = |, using the GA) 

per-user channel capacity curve close to the extreme upper 
bound, i.e., 1 bit per second per user The minor difference 
between the curve of ED and QD method justifies the accurate 

As is expected in 



approximation of QD in (lOi and (11 



Subsection |II-C the MD criteria capacity performance is close 
to ED and QD methods for high Eb/No values, which is shown 
in Fig. [4] Due to the fact that the ED method is closed to QD 
but with less complexity, we will choose the ED as the distance 
measure to compare with other criteria such as BER and sum 
capacity criteria. 

2) Comparison with Criteria: Fig. |5] demonstrates the 
per-user capacity of different proposed matrices optimized by 
BER, capacity, and ED criteria for different E^/Nq values. 
In addition, we compare these results with the Welch Bound 
Equality (WBE) codes which are introduced in pO) and pT) . 
Since WBE is optimum for Gaussian input distribution, there 
is no guarantee to be optimum for binary input vectors. Among 
the proposed matrices, the results of the ED method are close 
to that of BER scenario; this verifies that our approximation 
in (j9]) and ( [T2] i is accurate. As discussed earlier, to compute 
the BER criterion, a large array of bits needs to be processed 
which makes the computation of this criterion more complex 
compared to the ED criterion. Fig. |6] shows the per-user 
capacity of optimized matrices based on the capacity method, 
ED, and BER method versus loading factor (/3 = ^). in this 
figure we fixed an E^/Nq value (8dB) to compute the sum 
capacity. The performance of the capacity method decrease 
slower than other methods which shown the robustness of this 
method against loading factor. 



B. Results of GA for Binary Valued Signature Matrices 

Although, our proposed methods are applied to real valued 
matrices, we can use them to find sub optimal binary (±1) 
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GA). 



for various Ef,/No values (/3 = ^, using the PSO). 



matrices. Binary matrices are much simpler, in implementa- 
tion, than real valued ones. We consider an 4 x 5 binary matrix 
( at Table [Till ) which is optimized by the capacity criterion 
and compare it with another binary matrix derived from the 
ED criterion {A3 at Table IV 1. Fig. |7] shows per-user capacity 
of these two binary matrices with a real valued sub-optimum 



matrix derived from the sum capacity criterion (A4 at IIIi. 
This figure shows that the binary matrices developed by sum 
capacity criterion can be close to the real valued matrices 
however the binary matrix derived from the ED method is 
not as good. 



C. Results of PSO 

In this subsection, we present the simulation and numerical 
results of the PSO. Due to poor result of this algorithm as 
appos to GA for binary matrices, we only show the results for 
real valued signature matrices. 

Fig. |8] shows the performance of the sub-optimized matrices 
derived by PSO for the case when /3 = |. The curves of the 
BER, MD, and the capacity criteria are near the upper bound 
(1 bit per sec per user). In this figure, the results of the BER 
and the ED methods are very close to the capacity method 
unlike the GA depicted in Fig. [5] 

Fig. |9] shows the per user capacity versus the loading factor 
for various criteria. As expected, the performances of all 



the cases decrease with increasing f3. Similar to the GA the 
Capacity and BER criteria are the best criteria for large f3. 



D. Comparison of GA and PSO Results 

In this subsection, we evaluate the sensitivity of the pro- 
posed criteria with respect to the loading factor and opti- 
mization algorithms (GA and PSO). Fig. 10 demonstrates a 



comparison between the GA and the PSO for all the criteria 
(/3 = |). This figure shows per-user capacity of GA minus 
the PSO algorithm for various. Also, Fig {TT] presents the same 
results when /3 = |. A comparison between fig 10 and Fig. 
11 shows that for low loading factors (/3) and small values 
of Eb/NO the PSO performs better than the GA. On the the 
hand, for high values of betta the GA performs better Also, 
these two figure show that for the ED and BER criteria the 
choice of the GA and the PSO algorithms do not make any 
differences. 

Note as discussed in Section II the PSO algorithm is about 
5 times faster than the GA algorithm. Although note that, in 
general, the GA algorithm performs slightly better that PSO 
for real signature matrices but it is the only choice for binary 
matrices. 

In the following section we propose a method to design 
large size signature matrices for increasing the sum capacity. 
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Fig. 9. Per-user capacity of sub-optimum matrices versus loading factor Fig. 11. Per-user capacity derived from the GA minus that from the PSO 



(Et/No = 8dB, using the GA). 



for various measures versus the Ef, /Nq (f} = | ). 
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Fig. 10. Per-user capacity derived from the GA minus that from the PSO 
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V. Designing Large Signature Matrices 



The GA and the PSO algorithms are not suitable for large 
signature matrices such as 64 x 128. We thus propose a 
method of enlarging signature matrices from a sub-optimal 
small scale signature matrix. This section presents the method 
for enlargement from sub-optimal signature matrices. 

We propose to derive optimum matrices in a different way. 
Instead of the proposed methods discussed in Section II we 
derive sub-optimal signature matrices from small values of n 
and m and then by tensor products we construct a much larger 
signature matrix called an enlarged matrix. 

Instead of directly deriving optimum matrices using the 
proposed methods discussed in Section |ll] we derive an 
optimum signature matrix for a certain small value of n and 
m and later using "Kronecker product" yields sub-optimum 
matrix for the given dimension. We refer to this matrix as an 
"enlarged matrix". This section presents the theoretical proof 
for the practicality of the proposed enlarging method as well 
as the proposed ML decoder design. In this paper we manly 
consider signature matrices for binary and non binary cases 
[part II and Part III]. Here we try to extend these results real 
signature matrices for any types of inputs. 



A. Design Procedure 

The following theorem provides the necessary mathematical 
ground for enlarging signature matrices. 

Theorem 1: Suppose that A is a real mxn signature matrix 
for a CDMA system with binary or non binary input and 
G is a fc X fc reversible matrix such that its columns are 
normalized. Also, assume C{kn, km, (T]y\A) is the channel 
capacity assigned to matrix A in the presence of AWGN noise 
with variance ctjv. Denote (g) as the Kronecker product and let 
B = G (Ki A, then. 



(13) 



The equality holds if and only if G is unitary. 

Proof: Considering ([T]i as the model of CDMA systems 
and B as the signature matrix, we have: 



Ykmxl — ^Xknxl + Nkrnxl 



(14) 



where N is the noise vector with variance o-jy. For an arbitrary 
vector U, define Ui as the i*'* segment of U containing m 
entries; thus U'^ = [Ui,U2, . ■ . ,Uk]- Multiply both sides of 
( 14]) by D = G"^ (g) and define Z = BY, hence: 

Z = (G-i ® U)(G ® A)X + (G-i ® I^)iV 

From the Kronecker product properties, we have: 

(G-^ (g) I^)(G «) A) = G-^G (g) I„,A = Ik «) A 

Denote M = DiV. Thus for each i, the entries of Mi 
are independent Gaussian random vectors with variance 
<^M,i — \\G~'''{i,:)\\2(JN- Having h{Z), we can calculate 
h{Y) as follows: 



h{Y) = h{Z) - log2 |detD| 
From the Kronecker product properties: 

det (G-^ (g) Im) = (dctG-^)"(detI„,)'^ 
Considering ( [T5] l and ( [T6| , we have: 

h(Y) = + mlog2 |detG| 
Similarly, for N and M, we have: 

h{N) = h{M) + m log2 I det G| 



(15) 



(16) 



(17) 
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Bear in mind (|4]), we can write, 

C{kn,km,,aN\B) = max h{Z) - /i(Af) (18) 

P{X) 

Since Z — [Zi, Z-i^ - ■ ■ , Z^]'^ , we have the following upper 
bound for h{Z): 



3) Join X^'s to construct X such that X = 



\X'\X' , 



.,X>^']. 



The following is an example the Following example clarifies 
this algorithm: 



h{Z)<J2HZ^) 



(19) 



i=l 



III 



m 



Considering the relation between Y and Z and by using Q, 
we can derive the following equation: 

maxh{Zi) = C{n,m,aM,i\-^) + ^ log2 27recr^ - TOlog2 ||G' 
p{x) I 

where G^^(i, :) is the z*'* column of G^^. From Q for Y , 
we have: 

fc 

C(fcn,fcTO,CTjv|B) < ^C(n,TO,(||G-i(i,:)||(Tw)|A) 
i=i 

IdetGI 
+m logo — r 

ntil|G~Mv)ll 
If G is unitary, the equality holds and we have: 

C(fcn, /cm, CTArjB) = kC(n^m^a 

Assume, G is not unitary. Since columns of G are normal- 
ized, IdetGI < ntil|G(:,«)|| = 1 and \\G,-^{%,:)\\ > 
1. Thus, due to the reverse relationship between the sum- 
capacity and noise power, C(n,m, (||G^-'-(i, :)||(Tjv)|A) < 
C{n,m,aN)\A) and, 

C{kn,km,aN[B) < A:C(n, m, cr^vIA) 

Therefore, C(kn, km, CTAr|B) is maximized when G is unitary. 

■ 

According to this theorem the best choice of G for enlarging 
a signature matrix is an unitary matrix. In practice, we use a 
binary Hadamard matrix with normalized columns as G. Let 
Hji be a Hadamard matrix of order n, then H2 iS) Hn is a 
Hadamard matrix of order 2n. Using this fact, we can enlarge 
a signature matrix by a factor 2*^ where k = 1,2,... and 
construct Hadamard matrices from only 7?2- It is important 
to mention that although the differences between capacity of 
different criteria may not be noticeable, for large scale matrices 
it becomes important. 

The enlarged matrices can be decoded using a new simple 
ML decoder with significant reduction in the complexity of 
ML decoding. This decoder was originally introduced by Q 
for binary matrices and its none binary version discussed in 
Q and in this paper we generalize it to real matrices. 

B. The new ML decoder Algorithm 

The following is the step by step procedure of this decoder: 

1) Define Z = VkH,^^ ® I-mY and split it into Z^j^i' * = 
1,2,..., A:. 

2) Decode Z*'s using "tensor decoder" ||4j and obtain X^'s. 



Example 1: Suppose that the enlarged signature matrix is 

D = 4= Ho ® A4X5, where A4X5 is A4 from Table 

v 2 

the appendix; consequently, 

1 

71 



D 



Assume that X — 11,1 



A 
A 



A 

-A 



Y = 'DX + N 



1,-1,-1,-1,-1,1,1,-1]^, then 

-1.4586" 
-0.5227 
-0.8251 
-1.3148 
0.9584 
-0.1522 
3.7170 
2.0180 



where is an 8 x 1 AWGN vector. For decoding at the 
receiver, we split Y into two equal length vectors and 
y2 such that Y = [r^, F^jt ^^^^^ j^^^g^ 



y2 



[-1.4586, -0.5227, -0.8251, -1.3148]^ 
[0.9584, -0.1522, 3.7170, 2.0180]"^ 



Define: 



Z' - 
Z^ = 



y2] 



1 ryi _ y2] 



\/2 



The decoding of Z^ , Z^ using the tensor decoder yields X^ — 

[1. 1, -1, -1, -1]^ and X^ = [-1, -1, 1, 1, -1]^ where X = 

[x\xY. 

In general, for an km. x kn signature matrix which is 
enlarged from an m x n matrix, the usual ML decoding 
needs 2'^"'^''" Euclidean distance measurements while the 
new decoding needs fc2"^™ Euclidean distance measurements. 
As an example, suppose that we have a 64 x 80 signature 
matrix which is enlarged from a 4 x 5 signature matrix. The 
ML decoding for this matrix needs 2®^ * 2^*^ = 2^*^ Euclidean 
distance computations while using the "tensor ML decoder", 
we only need 16 * 2 = 32 Euclidean distance measurements. 

VI. Conclusion and Future Work 

In this paper, we have defined and derived the capacity for 
a specific signature matrix and a given Eb/No for overloaded 
CDMA systems. In order to find sub-optimum matrices, a 
number of optimization criteria have been introduced; namely, 
the capacity, BER and the distance criteria. We have modelled 
these criteria and shown that they differ with each other in 
terms of optimality and complexity. Our simulation results 
demonstrate that our proposed ED criterion presents the best 
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performance regarding the complexity and optimality accu- 
racy. To derive the sub-optimum signature matrices, we applied 
the GA as well as the PSO algorithms. In addition to the 
real valued matrices, we have also applied our methods in 
finding sub-optimum binary signature matrices. For large scale 
systems, instead of directly optimizing signature matrices, we 
enlarge small sub-optimal matrices using Kronecker products. 
We have shown that the capacity of these enlarged matrices is 
increased by the enlargement factor k. We can employ simple 
ML decoding for such enlarged matrices which significantly 
reduces the implementation complexity while maintaining 
optimality. While in this work, we have applied the GA and 
the PSO algorithms for optimization purposes,in order to avoid 
local minima, we suggest to work on "simulated anneaUng" 
as future works. Furthermore, we suggest derivation of sub- 
optimum matrices for non binary discrete valued signatures. 



[15] R.Hassan, B.Cohanim, O.de Week, G. Venter, "A Comparison of Particle 
Swarm Optimization and The Genetic Algorithm," 1st AIAA multidisci- 
plinary design optimization specialist Conf., No. AIAA-2005-1897, Austin, 
TX, 2005. 

[16] J. Brewer, "Kronecker Products and Matrix Calculus in System theory," 
IEEE trans, on Circuits and Systems, vol.25, no.9, pp. 772- 781, Sep 1978. 
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TABLE m 

Sub-optimum matrices based on the capacity method using GA,(note: A2 is optimized for Ei,/No 

OPTIMUM FOR Eh /No = SdB). 



lldB and other matrices are 



Capacity method 



A4 



0.1235 0.3177 0.7605 0.8739 0.4069 
0.3723 0.9240 0.5021 0.0154 -0.2553 
0.9764 0.3895 0.7448 -0.9375 
-1.0000 0.1711 0.4241 0.6451 
0.8529 0.6424 0.0930 1.0000 
■ 1 1 0.969 0.468 l' ' 

0.424 -1 0.5 -0.871 0.5 

1 0.015 -0.906 -0.75 0.719 
0.430 0.995 -0.938 0.984 0.984 



-0.3724 0.7299 
-0.5254 0.2538 



As 



-0.0115 
0.9584 



1.0000 0.5408 
-0.3117 -0.7224 



"1 


1 


1 


1 


1 " 


1 


-1 


1 


-1 


1 


1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


-1 



TABLE IV 

Sub-optimum matrices based on the ED method using GA,(the matrices are optimum for Ei,/Nq - 



ED method 



0.0591 0.8787 
-0.9198 0.1760 



-0.6226 0.4163 0.2166 
0.1907 0.6094 0.8851 



0.9572 0.4704 
-1.0000 0.8393 
0.3995 0.6776 



0.5922 
0.3621 
-0.7468 



0.1288 
0.7090 
-0.1777 



"-1 


-1 


1 


-1 




-1 


1 


-1 


1 


1 


1 


1 


-1 


-1 


1 


-1 


-1 


-1 


-1 


-1 



TABLE V 

Sub-optimum matrices based on the QD, MD, and BER methods using GA,(the matrices are optimum for Ei,/No = SdB). 



Size 



QD method 



MD method 



BER method 



0.4520 -0.3740 0.9029 0.1059 
-0.7780 0.3048 0.9585 -0.6561 
0.9163 0.4018 0.3265 -0.0717 
0.5315 0.9989 -0.9456 0.5273 0.4257 
0.4364 0.3203 0.5859 -0.9514 0.7039 



0.5924 0.1238 0.4630 -0.4371 
0.0557 0.4388 0.6436 0.5020 
0.9595 0.4137 0.0075 0.4325 



0.2502 0.4917 
0.6206 0.9009 
0.9903 0.2592 



0.1048 
-0.9958 
0.4383 



-0.9300 
0.4022 
0.9961 



TABLE VI 

Sub-optimum matrices based on various methods using PSO,(the matrices are optimum for Ei^/Nq = 8dB). 
Mthod 2 by 5 3 by 4 



Capacity 




'l.OOOO 1.0000 1.0000 -0.3120" 
0.9419 1.0000 -0.6067 0.0812 0.6859 




1.0000 0.3137 
1.0000 1.0000 
1.0000 1.0000 




ED 




'l.OOOO 1.0000 1.0000 0.5432 -0.0269" 
0.5206 0.1099 -0.2031 1.0000 1.0000 j 




1.0000 0.0665 1.0000 
1.0000 1.0000 
1.0000 1.0000 




MD 




0.3045 0.6719 1.0000 0.2925 -0.0804 
1.0000 0.2708 0.0711 -0.7045 1.0000 




1.0000 1.0000 0.0483 
1.0000 1.0000 0.0574 
1.0000 0.0701 1.0000 




BER 




'1.0000 -0.7644 1.0000 0.4113" 
1.0000 1.0000 0.5402 -0.4707 1.0000 




1.0000 1.0000 -0.2516 -0.9898 
1.0000 -0.1536 1.0000 -0.0209 
1.0000 -0.6976 1.0000 





